Foundation of Fractional Langevin Equation: Harmonization of a Many Body Problem 
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In this study we derive a single-particle equation of motion, from first-principles, starting out 
with a microscopic description of a tracer particle in a one-dimensional many-particle system with 
a general two-body interaction potential. Using a new harmonization technique, we show that the 
resulting dynamical equation belongs to the class of fractional Langevin equations, a stochastic 
framework which has been proposed in a large body of works as a means of describing anomalous 
dynamics. Our work sheds light on the fundamental assumptions of these phenomenological models. 
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I. INTRODUCTION 

The stochastic dynamics of many-body systems with 
general two-body interactions are inherently difficult to 
solve. There are, however, a few idealized exactly solv- 
able models that have served as benchmark cases from 
which collective effects have been deduced [TH3]. One 
example is the one-dimensional over-damped motion of 
non-passing hard spheres (so called single-file diffusion) 
in which a tracer particle behaves sub diffusively [4] [5]. 
Processes displaying anomalous diffusion occurs in a 
range of many-body systems [6-8 , especially in biology 
[9j[10], and are, apart from a few exceptional cases, mod- 
eled in phenomenological ways. 

The one-dimensional motion of identical Brownian par- 
ticles (BPs) which are unable to pass each other is well 
studied theoretically fTTHT5] . A tracer particle in such 
single-file system exhibits sub-diffusion; its mean square 
displacement (MSD) is proportional to t 1 ! 2 indicating 
slow dynamics [11] [16]. There exists a wide range of 
experiments on tracer particle dynamics in diverse sys- 
tems which show the t^-behavior. Examples include 
colloids in one-dimensional channels [17-19 , an NMR ex- 
periment involving Xenon in microporous materials [20] , 
molecular diffusion in Zeolites [2TJ [22] , moisture expan- 
sion in ceramic materials [23] , and a study of Ethane in 
a molecular sieve [24] . 

Independent of the developments in the field of single- 
file diffusion, recently a stochastic fractional Langevin 
equation (FLE) has gained much interest [25-27 . In the 
FLE a derivative of fractional order a replaces the usual 
first order time derivative in the overdamped Langevin 
equation (d/dt — » d a /dt a ). FLEs with 0<a<l are able 
to describe phenomenologically a range of physical phe- 
nomena such as motion in viscoelastic media [28] . In the 
presence of a binding harmonic field, Xie and co-authors 
used the FLE to model protein dynamics [26j [2£ and 
a = 0.51±0.07 was deduced from experimental observa- 
tions. Here we derive an FLE with a = 1/2, in the pres- 
ence of an external force field, starting from a many-body 
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FIG. 1: Illustration of the many-particle system studied here, 
(a) Particles interact via the two-body potential V(\x n (t) — 
x n '(t)\) and cannot pass each other, (b) The harmonization 
procedure amounts to mapping the interacting particle system 
onto a harmonic chain in which the interactions are captured 
by the effective spring constant k. 



theory. Thus we show that the 1/2-FLE is expected to 
be universal for a large number of experiments describing 
interacting tracer particle dynamics. 

Usually single-file models consider BPs with hard-core 
interaction. These models can be mapped onto a non- 
interacting system with known methods p~2j [30] . In ex- 
periments the interaction of BPs is hardly ever hard-core. 
In this paper we consider rather general two-body inter- 
actions between BPs. We present a new method to deal 
with this many-body problem, which we call harmoniza- 
tion. With this method we are able to effortlessly derive 
many previous results, for example the Kollmann rela- 
tion for the MSD [14 , to justify from first-principles the 
FLE, and to derive many new results, such as the dis- 
tance correlation function between a pair of particles. 



II. SYSTEM DESCRIPTION 

We consider particles undergoing one-dimensional 
over-damped Brownian motion in an infinite system 
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where two particles n and n' interact via the two-body 
potential V(\x n (t) — x n f(i)\), where x n (t) is the position 
of the nth particle. The potential has a hard-core part, 
so the particles cannot pass each other, but otherwise 
it is rather general [31]. The Langevin equation for the 
motion of particle n reads 

= £ f M*) " Xn ' + ^(V + fn(t), (1) 
n' 

where £ is the friction constant (D = k B T/^ is the free- 
particle diffusion coefficient), f = —dV/dx n is the force 
due to interactions with the surrounding particles, r] n is 
a zero mean white Gaussian noise with (rj n (t)rj n f (t')) = 
2^k B TS(t — t')5 nn > where k B and T are the Boltzmann 
constant and temperature. f n is an external force. The 
particles are initially taken to be in thermal equilibrium. 
Our main interests are in the dynamics of a tracer par- 
ticle position and in distance fluctuations between parti- 
cles. These quantities are, however, intractable from the 
many-body problem for a general V. Therefore, we 
introduce a new technique - harmonization. 



III. HARMONIZATION TECHNIQUE 

The philosophy behind harmonization is to map the 
original system (A), i.e., a system described by Eq. ([I]), 
to a system (B) consisting of beads interconnected by 
harmonic springs. Consider two particles in system A 
with N ^> 1 particles in between, and which, in equi- 
librium, are at an average distance L eq from one an- 
other. We let T{L, N) be the extensive free energy due 
to the N particles at a fixed temperature T. Fluctu- 
ations of L around L eq are small. Hence we may ex- 
pand F(L,N) ~ F(L ec{ ,N) + k N (L - L eq ) 2 /2 where 
we introduced a macroscopic spring constant: k^ = 
d 2 J r (N 1 L)/dL 2 \L eci . Now we replace system A, by a sys- 
tem B of beads connected to their two nearest neighbours 
by springs with spring constant n. We relate the spring 
constants n and k^ by requiring that the total free ener- 
gies associated with identical displacements from L eq in 
the two systems are the same. This gives 
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d 2 J r (N, L) 



dL 2 



(2) 



To see this note that for N particles in system B the to- 
tal free energy change is Nn(L — L eq ) 2 / (2N 2 ) [since the 
displacement of one spring is (L — L eq )/N and we have 
N such springs] which is equal to the free energy of a 
single macroscopic spring in system A, kw(L — L eq ) 2 /2, 
when Eq. ([2| holds. The effective spring constant n is 
an intensive thermodynamic quantity, which can be ob- 
tained from the original system's equation of state. From 
the definition of pressure P = —(dF/dL)^ and isother- 
mal compressibility xt — — (dL/dP)N/L, one obtains 
k = -N(dP/dL) N = p/xr where p = N/L (AT, L -> oo) 



is the particle number density. For a one-dimensional 
gas of hard-core interacting point-particles the equa- 
tion of state is given by P = Nk B T/L which leads to 
k = p 2 k B T. Similarly, for systems consisting of 6-sized 
hard rods (Tonks gas) we have from the van der Waals 
equation P = Nk B T/(L-Nb) that n = p 2 k B T(l- pb)~ 2 . 

The final step in the harmonization method is to re- 
place the non-linear two-body interaction in Eq. by 
forces from the nearest-neighbor spring coupling, i.e., 

dx (t) 

(3) 

We above implicitly assumed that there has been time for 
particles to interact with neighbors, t ^> T[ nt — l/(p 2 D). 
Equation ([3| will be justified later, via simulations, and 
below by showing agreement with known results for 
single-file systems. Under the assumption t ^> T{ nt we 
can take the continuum limit and turn x n (t) into a field 

„dx(n,t) d 2 x(n,i) 



+ /(n,t)+*y(M). (4) 



dt dn 2 
This relationship is our harmonization equation from 
which previous exact results are recovered and new ones 
derived. In the following subsection the MSD of a tracer 
particle is discussed. We note that a mapping similar 
to Harmonization was applied to the simple exclusion 
process in [5]. Here we consider a general two body in- 
teraction showing precisely how to compute the effective 
spring constant from equilibrium concepts (i.e., the com- 
pressibility) . 



A. Mean square displacement 

Consider the case of no external forces, f(n,t) = 0. 
Equation (ij) is then equivalent to the Rouse model from 
polymer physics (see for instance [32 J. In the following 
we consider a tracer particle labeled T and consider its 
MSD. We will arbitrarily choose the particle n = with 
position x r (t) = x(n = 0,t). The MSD, {5x 2 r (t)) = 
{[xq-(t) — xj-(0)] 2 ) (() denotes average over noise and 
initial conditions), is [33] 



(Sx 2 r (t)) = k B Ty/At/(^K). (5) 

The derivation is relegated to Appendix [Aj If k for the 
gas of 6-sized hard rods is used one obtains (5xj*(i)) = 
p^il-p^^Dt/^, which agrees with gj HSJ [34]^] . In 
Ref. [14 Kollmann showed that (Sxj-(i)) ~ t 1 ! 2 regard- 
less of the nature of interactions as long as mutual pas- 
sage is excluded. In p articular Ko llmann derived the re- 
lation (5xj-(t)) ~ So^4:D c t/(p 2 7r), where 5o is the static 
structure factor at zero wavevector and D c is the collec- 
tive diffusion constant. Equation (5]) gives via Eq. ^ 
a relation between the MSD and the free energy of the 
system, while the Kollmann relation relates the MSD to 
physical observables 5o and D c . Equivalence between 
our results and [14] is found using 5q = k B Tpxr and the 
relation D c = 1/(pxtQ [37] - 
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IV. FRACTIONAL LANGEVIN EQUATION 1Q 2 



Using our harmonization technique we recovered 
known single- file results. Now we take the harmoniza- 
tion one step further and derive an FLE for the posi- 
tion of the tracer particle xj-(i) = x(n = 0,£). Tak- 
ing the Fourier and Laplace transforms [38] x(q, s) = 
dn J °° dt e- iqn - st x(n, t) of Eq. (E| gives 



x{q,s) 



ri(q,s)+Zx(q,t = 0) + f(q,s) 
£s + nq 2 



(6) 



Note that we for all functions indicate a Fourier trans- 
form with the variable q and a Laplace transform by the 
variable s. Subtracting 27r5(q)x<r(t = 0)/s from both 
sides of (|6|, rearranging, and taking the inverse Fourier 
transform at n = yields 

7(«)W«) " Mt = 0)] = f)f{s) + f T (s), (7) 



where 7(5) = ^/A^k/s is a fractional friction ker- 
nel, and the bar over a quantity means yr( s ) — 
dn exp( — y/^s/~K,\n\)y(n, s). The effective noise is de- 
fined as 7? eff (n, s) = 77(71, s) + £[x(n, t = 0) — xj-(t = 0)] 
and includes Gaussian noise as well as randomness in the 
initial conditions relative to the tracer particle. If the 
external force acts only on the tracer particle, /(n, t) = 
fj-(t)5(n), the inverse Laplace transform of Eq. (TtI) yields 



-d x ' 2 x r (t) 
1 dt 1 / 2 



(8) 



where we introduced the Caputo fractional derivative 



d a f(t) 
dt a 



1 



T(l - a) J \t - t 



dt' df(t') 



dt' 



(9) 



of order a = 1/2. Equation ([8| is the sought for FLE 
with fractional kernel j(t) = ^/4£ft/ (irt) [39] and is in 
agreement with the long time limit of the result proposed 
phenomenologically for hard-core interacting point parti- 
cles in [40]. Notice that the form j(t) oc t~ x l 2 is a direct 
consequence of the harmonic expansion and the assump- 
tion of over-damped dynamics. Assuming thermal ini- 
tial conditions one straightforwardly (see Appendix [B| 
shows that the effective noise satisfies the fluctuation- 
dissipation relation 



(fjf(t)fjf(t')) = k B T 7 (\t - t'\). 



(10) 



This was expected since the Langevin equation was orga- 
nized with the external force conjugate to xq-(t) as a term 
on the right hand side (see [41]). For a constant force, 
fr(t) = F, one can deduce from Eq. ([8| (see Appendix 
[Cl the generalized Einstein relation 



(x r (t)) F = F(Sx 2 T (t))/(2k B T), 



(11) 



where {xj-(t))p is the average shift in position in the pres- 
ence of the force and (Sxj-(i)) is the MSD in the absence 




FIG. 2: Autocorrelation function obtained from molecular dy- 
namics simulations of 10 4 point particles on a ring (solid lines) 
with corresponding theoretical predictions (dashed lines). 
Upper black lines: tracer particle in an external harmonic 
potential with the Mittag-Lemer function ( 13 ) (parameters 
are: p 0.25, k B T = 1, f = 0.5, moo 2 0.0158, resulting 
in r ~ 501). Lower red lines: inter-particle distance correla- 
tion function with the prediction of Eq. ( 18 ) (parameters are: 
p = 0.4, k B T = 1, £ = 0.5, ria - n d = 10, giving r d i s t ~ 313). 



of the force, i.e., as given by Eq. ([5|; this Einstein rela- 
tion generalizes the results in [42] to systems with general 
V. For the case of a periodic force fr(f) — i^ocos(cjo^), 
we find asymptotically at long times a non-trivial 45° 
phase-shift between the applied force and mean displace- 
ment (see Appendix O) : 



(x r (t)) Fo 



F 



cos (ujot — 7r/4) . 



(12) 



The response of a tagged single-file particle to a harmon- 
ically oscillating force was previously obtained in a differ- 
ent way in [40] . More generally it has been obtained from 
the starting point of the fractional Langevin equation in 



A. Tracer particle in a harmonic potential and 
simulations 



One of the predictions of FLE theory is that the au- 
tocorrelation function (x j-(t)x j-(0)) , under the influence 
of an external harmonic potential, decays as a Mittag- 
Lemer function [25j [26] . For the tracer particle in ther- 
mal equilibrium with respect to the harmonic force fj- = 
—muj 2 xq-(t) we find from Eq. ^ 



(x r (t)x r (0)) 

where E 1/2 (-(t/r)^ 2 ) 
Leffler function and r 



£i/ 2 (-(t/r) 1/2 ), (13) 



e^erfc^/r) 1 / 2 ] is the Mittag- 
= A^n/(muj 2 ) 2 . For t > r this 
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leads to a decay (xj-(t)xj-(0)) ~ k B T \Jr / (nt) / \muj 2 ) . 
We tested the prediction of Eq. ( p~3| ) numerically by sim- 
ulations of hard-core interacting point particles. The re- 
sult is shown in Fig. |2j agreement with the analytic 
prediction is excellent (without any fitting) in the time 
regime where particles have collided and harmonization is 
valid, i.e., when t ^> ri nt = l/(p 2 D). The slight deviation 
at shorter times is in accordance with the interaction-free 
Brownian motion of the simulated particles prior to col- 
lisions, see also [40]. An autocorrelation function with 
Mittag-Lemer decay of index 1/2 was recorded in the 
experiments [29]. The harmonization procedure can be 
applied to other problems than tracer particles. This is 
illustrated below. 



V. INTER-PARTICLE DISTANCE 
CORRELATIONS 

Donor-acceptor data from conformational dynamics of 
proteins [44] was recently modeled using an FLE with a 
harmonic potential [i.e., Eq. ([8| with fj- harmonic] [26] . 
Fracton models [45] and FLEs based on the Kac-Zwanzig 
model [46] have also recently been studied. Contrasting 
more phenomenologically oriented approaches, our har- 
monization technique allows us to attack problems re- 
lated to inter-particle dynamics on a first-principle level. 
In fact, considering inter-particle distance dynamics we 
will now show that: (i) a harmonic potential arises natu- 
rally and is not due to an external field as assumed so far, 
and (ii) the governing equation is a generalized Langevin 
equation (GLE) with a power law memory kernel which 
leads to anomalous relaxation, rather than an FLE. 

Defining x d [ st (t) = x(n a ,t) — x(n dl t) as the distance 
between an "acceptor" particle n a and "donor" n dl an 
equation for Xdist(^) is obtained by subtracting Eq. ^ 
for x(n dl t) from the corresponding one for x(n a ,t). If 
external forces acting only on particles n a and n d are 
considered, f(n,t) = f na (t)S(n - n a ) + f nd (t)S(n - n d ), 
we find 

7(5) [sadist (s) - x dist (t = 0)] = rf ff (n a , s) - ?f ff (n d , s) 

+ (/n.00 - /nato) (l " e~^ k ^) , (14) 



in Laplace space where k& 



Eq. (14) by 2(1 - e 



ist 



/sK,£/k d 



ft/l^a — n d|- Dividing 
ist ) will result in the last 
term on the right hand side becoming a force f d [ st = 
(/n a — /n d )/2 which is conjugate to the coordinate x d i st . 
Thus we rearrange Eq. ( [14] ) to the form of the GLE: 

dt' K{t - i>) dxA ^P = mist (t) - U'(x dist ) + f dist (t), 

(15) 

with friction kernel and a noise in Laplace space given by 

l(s) fcdist 



/ 

Jo 



JC(s) = 

^dist(s) = 



2(1 — e -v^<Adist) s ' 

fj eff (n a ,s) - rf ff (n d ,s) _ fedist 
2(1 — e-Vswf/fedist) s 



(16) 



Sx d i St (t = 0), 



where Sx d [ st (t) = x d [ st (t) — (n a — n d )/p. The constant 
terms subtracted in the definitions of the friction and 
noise are combined in the harmonic force — W(#dist) = 
-fcdistfcdistW - (n a -n d )/p\. The spring constant fe dist 
corresponds exactly to the |n a — n d \ ft- strength springs 
that are connected in series in between the donor- 
acceptor particles after harmonization. A straightfor- 
ward but lengthy calculation (see Appendix [D| confirms 
that (7?dist(^dist(0) = k B TJC(\t - t'\), as required by 
the fluctuation-dissipation theorem. One can invert the 
kernel /C(s) exactly and express it as a Jacobi theta func- 
tion 



K(t) 
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1 + 



(o,« 



-TdiBt/4* 



k, 



dist i 



(17) 



where #3(2, u) = 1 + 2 J2m=i uW cos (^ mz ) an d r dist = 
^Adist- Note tnat r dist = (\n a - n d \/p) 2 /D c so that 
fdist can be interpreted as the time it takes for the in- 
formation about the motion of particle n a to diffuse 
to nd- Examining Eq. (17) for t <C Tdi s t one finds 
JC(t) ~ (irk) which implies that the two particles 

do not influence e ach other . For longer times t ^> Tdist 
we have JC(t) ~ yj ^/(Ant). The factor 1/2 difference in 
the prefactor of JC(t) at long and short times means that 
the autocorrelation function will not decay exactly as a 
Mittag-Lemer function. Instead from Eq. ( 15 ) we obtain 



(<fadist(*)<fadist(0)) 



k B T 

^dist 



erf 



Tdist 

4t 



1 



t/4t 



V^r dist /4t 



(18) 

with (Sx dist (t)Sx dist (0)) - k B Tk,l t Jr d[s t/(47vt) at long 



times. In Fig. [2] we compare Eq. (18) to simulations and 
find excellent agreement (without fitting). 



VI. SUMMARY AND CONCLUDING 
REMARKS 

Throughout this paper we have shown that our harmo- 
nization technique can reproduce known results as well 
as providing new ones. But how come it works so well? 
When equation Q was obtained, a quasi-static approx- 
imation was used in the sense that the effective spring 
constants were calculated based on the equilibrium prop- 
erties of the system. To see why this is physically reason- 
able one can argue as follows. The MSD of a tracer parti- 
cle in a single-file system is proportional to t 1 / 2 and it will 
therefore cross a system of length L in a time on the or- 
der of L 4 . This is considerably slower than the relaxation 
time of the whole system which scales as L 2 /D c . Thus, 
in the long-time limit a tracer particle only sees parti- 
cles which have had sufficient time to reach local thermal 
equilibrium [47] . This is the reason why equilibrium con- 
cepts like free energy work so well here. This implies that 
it is the one-dimensional topology and the single-file con- 
dition that leads to slow dynamics of a tracer particle and 
the possibility to map it onto a harmonic chain. The FLE 
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with exponent 1/2 is therefore expected to be found in 
a vast number of over-damped systems. Our framework 
can, however, be applied to particle motion in higher di- 
mensions. For instance, particles embedded in networks 
in which ordering is maintained. 

In summary, we have presented a harmonization tech- 
nique which maps a stochastic many-particle system with 
general two-body potentials onto a system of intercon- 
nected springs. The interaction potential was reduced 
through equilibrium considerations to only one param- 
eter: the spring constant n related to the compressibil- 
ity of the particle system [48 . We derived, from first- 
principles, an FLE which predicted subdiffusive (slow) 
dynamics of a tracer particle. Derived expressions agree 
perfectly with rigorous well-known results when they are 
available. Under the influence of an external harmonic 
force, Mittag-Lemer relaxation was found which was cor- 
roborated by simulations of a hard-core system. It would 
be interesting to test the harmonization technique further 
with simulations beyond this hard sphere model. The 
dynamics of inter-particle distance was also addressed 
and the harmonization approach predicted a GLE rather 
than, as previously suggested, an FLE. Unlike the or- 
dinary Langevin equation, which describes a Markovian 
process and which is usually derived for a massive par- 
ticle colliding with independent gas particles, the FLE 
exhibits long memories which, as we showed here, are 
due to the many-body nature of the underlying dynam- 
ics. Thus, fractional calculus enters through many-body 
effects which might be the reason why it took so long for 
a natural microscopic origin to be uncovered. 
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Appendix A: Mean square displacement of a tracer 
particle 



In this appendix we calculate the MSD of the tracer 
particle position xj-(t) = x(n = 0,t) defined as 



(Sx 2 T (t)) = ([x r (t) - x T (t = 0)] 2 ), 



(Al) 



under the assumption of no external force f(n,t) = 0. 
In order to find the MSD we calculate the correlation 
(x(q, s)x(q\ s')) which we find by multiplying Eq. ([6| by 



itself and average over the noise (...): 

{x(q, s)x(q', s) = A init (q, q' , s, s') + A noise (q, q' , s, s'), 



A init (q,q f ,s,s f ) = 

A n0 ise{q^ q i s i s ) 



is 2 (x(q,t = 0)x(q',t = 0)) 
(si/ + g 2 )(sV + g' 2 ) : 
= 1 (ri(q,s)ri(q',s')) 
k 2 (sv + q 2 )(s'v + q' 2 ) ' 



(A2) 
(A3) 



where v — £/k and the variable name q' (s f ) like q (s) 
implies that the corresponding variable have been Fourier 
(Laplace) transformed. Above we used that the initial 
positions are independent of the future thermal noise. 

Starting with A no i se we first note that the Fourier and 
Laplace transform of the noise autocorrelation function 
is 



(r](q,s)ri(q\s')) 
Thus we can write 



4*rtk B T5{q + q f ) 



(A4) 



^- 1 k B T5{q + q f ) 



, (A5) 



^-noisefe*? •> 5 ? s ) — / . ,w . ? / \/ / , 9 / \ ' 

(s + s')(s + q z /iy){s / + q l jv) 
Taking the inverse Laplace transforms of this we find 

e -q 2 \t-t'\/v _ e -q 2 {t+t')/v 



A no \ se (q, q', t, t r ) = 2n8{q + q' 



Kq 2 /(k B T) 



(A6) 

If we take the inverse Fourier transforms of the above 
equation and evaluate it n = n r and t = t' we find 



= k B T 



dq 1 



-2q 2 t/v 



2tt 



2t 



(A7) 



where we used dz z~ 2 [l — e~ az ] = 

We now proceed to evaluate Ai n it ■ Using the inverse 
Fourier transform J r ~ 1 {2a/(a 2 + ^ 2 )} = e _a l n las well as 
the convolution theorem we can write 



t (n,n / ,8,8 / ) 



du 



p —\n—u\y/su 

du' (A8) 



g— |n 7 — u \yj s' v 

x , (x(u, t = 0)x(u', t = 0)) 



2v^7 



V 



In the harmonic chain, the particles are in thermal equi- 
librium with respect to the potential 



(A9) 



from which the equilibrium density is Pequiiib. = 
e -v/{k B T)/z where Z = J[n m dx m ] e - u ^ kBT \ Thus, the 
particles' initial positions are Gaussian variables which 
we can express as 

Er=l A - n > 

x n (t = 0) = { n = (A10) 

-Ei n iiA_ r n<0 
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where we have chosen the coordinates such that xo(t = 
0) = 0. The expected values of the A r are 



(A r -p- 1 ) = 0, 
((A r -p- 1 )(A r , 



k B T 



(AH) 

S r y , (A12) 



where d r y is the Kronecker delta. From this we find 
{[x(u, t = 0) - x(0, t = 0)] [x(u', t = 0)- x(0, t = 0)]) 

7 rp 

= uu'p- 2 + min(|w|, |u'|)0(W) — , (A13) 



where is the Heaviside step function and min(a, b) 
is the smallest value of a and 6. Inserting this initial 
distribution in Eq. ( A8 ) and setting n = n' = leads to 



A- m i t (n = 0,n' = 0,s,s' 



^1 duf 

K JO JO 



du r 2 min(u, u') 



2^/VJv 



u 2 k B T 



2k 
v 2 k B T 



du I duu 
o ./o 



/"OO /-u 
/ du' / 

Jo Jo 



dui 



2K Jq Jq yj SV 

Using the inverse Laplace transform 

f p — Uy/sV I p — U 2 v/4:t 

C 1 ' ' 



, (A14) 



(A15) 



evaluating the second integral in each of the terms above, 
and setting t = t' gives 



Anit(n = 0,71 = 0,t,t) 



2 fc B T 



7T ft 



-2u 



= k B T J -^(2 -y/2), 

V 



(A16) 



wher e v = £/k was used. If Ai n i t is combined with Eq. 
(A7) we find the desired result for the MSD for a tracer 



particle in thermal equilibrium 

(Sxj-(t)) = Ai n i t (n = 0,n = 0,t,t) + A noi8e (n,n,t,t) 



k B Ty/At/(^K), 



(A17) 



which is the result mentioned in the main text. We note 
that if the particles initially had been placed equidis- 
tantly, x(n, t = 0) = n/p, with no randomness in the po- 
sitioning, then A{ n [t would have vanished and the MSD 
would have been smaller by a factor y/2. 



Appendix B: Fluctuation-dissipation relation for ?7 eff 

Here we will find the noise-correlation in Laplace space 
for the more general case (f/ eff (n, s)fj eS (n', s')) from 



which the fluctuation-dissipation relation follows as the 
special case n = n' . The more general case will be needed 
in Appendix [D] 

We will divide the noise into two parts 

ff*(n,8)=fj(n,8)+f?* t (n,8), 



r 

77(71, s) = / 

j — c 



(Bl) 

where the first part is related to the original thermal noise 

due-^'-'^u.s), (B2) 
and the second part is related to the initial positions 

/oo 
due -vMn-«| 
"OO 

x£[x(u,t = 0)-x(n,t = 0)]. (B3) 
For the first part we have the correlation function 

due-^ n ~ ul 



(rj(n,s)rj(n',s')) 

/oo 
du' 
-00 



e -V7Z\W-u'\ ^ML S (u - u'). (B4) 



Doing the integrals one finds 

(fj(n, s)fj(n' \s')) 



(s + s')(s - s') 



( 



Ve -^|n-n'|j ^ 



For the part of the noise correlation that comes from the 
random initial condition we have 



(77 init (n,s)7/ init (n / ,s / )) = £ 2 / due-^ n ~ u \ 

J —00 

/oo 
du' e -V^|n'-u'| fafa t = 0) _ a .( n> t = 0)] 
-00 

x [x(u\ t = 0)- x(n', t = 0)]) (B6) 
Using Eq. ( |A13[ ) one finds after a bit of calculation that 



2k B Tj 
s — s' 
1 



1 



I SV 
fsv\n—n' I 



I S'V 



/ s'v\n— n' I 



• (B7) 



Combining the two noise parts leads to the general for- 
mula 



2k B Ty/jk ( e -^\ n - n >\ J_ ( 

s + s' \y/s vs 7 ' 



fsv\n— n I 



- ( 7 (s) e -V^I"-n'l + 7 ( s ') e -V^\n-n'\^ (gg) 



S + S 

Setting n = n' finally gives 



s cff, >u , cff , ,i\ \ _ fc B r[ 7 (s) + 7 (a / )] 

s + s' 



(7? eff (n,s)7? eff (n,s')) = 



(B9) 



which is the Laplace transform of the sought for 
fluctuation-dissipation relation (ff^{n,t)ff^{n,t')) = 
k B Ti(\t-t\). 
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Appendix C: External force on a tracer particle 

Here we consider a force acting only on particle 0: 

f(n,t) = F(t)5(n). (CI) 

Taking the average of Eq. (|6| with respect to the 
zero-mean noise and using the inverse Fourier transform 
F~ 1 {2a/(a 2 + q 2 ) } — e _a l n l with the explicit expression 
for the force Eq. (CI) leads to 



(x(n,s)) = 



j — c 



du (x(u, t = 0)) 



e -yfsv\n-u\ F(s) e -yfsv\n\ 

2^J~sv 2k 



np F(s) e~^\ n \ 
s 2k \fsv 



(C2) 



Here we used that (x(n,t = 0)) = n/p. 



1. Periodic force /(n, t) = 5(n)F cos(cjo£) 

If we use the complex representation of the force 
f(n,t) = S(n)^t[Foe tUJot ] (where 5ft represents the real 
part) we have the Laplace transform 



f(n,s) = 5(n)$t 



F 



S — ZUJq 



(C3) 



Using this in Eq. ( C2 ) gives 



=n/p 



(x(n, s) — x(n, t = 0)) = 5ft 



-yfsv\n\ 



2K y/sv(s — IUJq) 



(C4) 



For n = 0, i.e we track the tagged particle on which the 
force acts, we have 



(x r (t)} = 5ft 



F e iuJot Eif(^u^i) 



2k 



(C5) 



where we used £ 1 {(y/s(s — a)) 1 } = e at Erf(\/ai)/y/a. 
Making the replacement v = £/k and i 1 ^ 2 = e* 71 "/ 4 leads 
to 



(x r (t)) = 5R 



e i(»ot-*/*) E rf( y /fa^ 



(C6) 



For large f we have Erf (\fTujQt) ~ 1 and therefore 



(* r (t)> 



cos (cjo^ ~ 



(C7) 



after taking the real part. 

2. Constant force f(n,t) = S(n)F 

The result for a constant force is obtained in the cjq 



limit of Eq. ( C6 ) . Using the short time expansion 



Erf(v / i^)~2\/^=2 A /^/ 4 

7T V 7T 



combined with Eq. ( C6 ) gives 



(x r (t)) = F 



TTK^ 



(C8) 



(C9) 



Together with Eq. ([5| this demonstrates the generalized 
Einstein relation, Eq. ( pTj ) in the main text. 



Appendix D: Fluctuation-dissipation relation for ?7dist 

Here we address the fluctuation-dissipation relation for 
77dist • Most of the work for deriv ing this was done in 
Appendix pi when deriving Eq. (B8). What remains 



is to work out the correlations of the part ^st( s ) 
&dist&£dist(£ = These turn out to be: 



kBTk dist 



k B Tk d 



^dist / 
' V 



1-e 



fsv\n^—ri a \ 



f s (n d ,s)^(s')) = -(fr(n„s)rffi£.(s')) 



Combining these correlations with Eq. (|B8j) one arrives 
at 

(%ist(s)r?dist(s')) = + > ( D1 ) 

which is the Laplace transform of the sought for relation: 
(Vdist(t)r ]dist (t , )) = k B TIC(\t-t , \). 
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